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1.0  INTRODUCTION 

Among  Ihc  requirements  encountered  in  ground-based  simulation  test  facilities  is  the 
need  to  determine  the  gaseous  composition,  temperatures,  and  How  velocity  in  both 
simulation  flow  fields  and  the  exhaust  plumes  of  various  engines.  A  variety  of  techniques  are 
being  developed  and  applied  to  such  problems,  but  the  measurement  of  the  properties  of 
minor  species  has  been  most  difficult.  This  capability  is  desired  to  provide  both  spatially  and 
temporally  resolved  measurements  of  the  number  densities  of  various  low  concentration,  or 
minor,  species  in  combusting  flow  fields  and  engine  exhaust,  plumes,  with  OH.  H,  and,  in 
some  cases,  H2  being  exemplary  species. 

It  is  known  that  for  very  low  concentrations  (mole  fractions  on  the  order  of  10~ 2  or  less) 
of  such  species,  local  measurements  using  Raman  scattering  are  impractical  because  of  low 
signal  levels.  However,  fluorescence  techniques  do  possess  the  necessary  requisites  of 
molecular  specificity,  and  spatial  and  temporal  resolution.  Also,  since  fluorescence 
excitation  cross  sections  exceed  Raman  scattering  cross  sections  by  over  eight  orders  of 
magnitude,  the  sensitivity  of  fluorescence  is  exceptionally  good.  Further,  unlike  charged - 
particle  fluorescence  excitation,  the  laser-induced  fluorescence  (I.IF)  technique  is  energy 
eigenstate  specific,  which  simplifies  considerably  the  resulting  fluorescence  spectra.  These 
advantages  of  LIF  arc  lessened  considerably,  however,  by  leeognilion  of  the  existence  oT 
quenching  processes  of  the  excited,  radiating  species,  for  now  the  radiation  intensity  is  no 
longer  a  linear  function  of  the  local  density.  Furthermore,  the  degree  of  nonlinearity  is  a 
function  of  the  number  density  of  each  species  in  the  sample. 

However,  it  is  known  that  for  sufficiently  intense  laser  beam  sources  interacting  with 
simple  two-  and  three-energy  level  systems,  the  laser -induced  stimulated  emission 
predominates  over  the  loss  ol  excited  state  species  by  collisional  quenching.  Consequently, 
for  such  saturation  conditions,  collisional  quenching  is  unimportant,  and  the  laser-induced 
fluorescence  provides  a  direct  measurement  of  the  local  species  density.  Unfortimally,  most 
molecular  species  of  interest  are  not  accurately  described  by  .simple  two-  and  iltrec-levd 
systems  but  rather  arc  multi-state  systems.  To  evaluate  the  feasibility  and  advantages  of  the 
use  ol  laser-induced  fluorescence  lor  such  species,  it  is  necessary  to  describe  analytically  the 
response  of  such  systems  to  intense  laser  beam  sources  and  to  determine  the  effects  and 
importance  of  collisional  quenching  processes.  This  report  describes  the  lesults  of  a  study  of 
the  time-dependent  fluorescent  signal  of  a  commonly  encountered  plume  gas  species  CO. 

During  this  sludy  it  became  clear  that  even  if  saturation  conditions  were  not  achieved,  a 
second  method  of  application  of  LIF  was  possible  if  the  fluorescent  species  was  a  trace  gas 
species  in  a  binary  mixture.  For  a  situation  in  which  collisional  quenching  of  the 
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fluorescence  ol  ihe  trace  species  exists,  the  fluorescence  signal  is  a  direct  function  of  I  lie  local 
density  of  the  host  gas,  or  major  species,  of  the  mixture;  i.e.,  for  a  specified  excitation 
source  the  fluorescence  response  of  the  trace  species  is  dominated  and  determined  by  the 
density  ol  the  quenching  species,  and  this  dependence  should  yield  a  direct  measure  of  the 
major  species  density.  To  investigate  the  possibility  of  measurement  of  the  density  of  the 
major  species  using  tracer  gas  impurity  fluorescence,  it  was  decided  to  focus  attention  on  a 
problem  area  for  which  such  a  technique  would  be  useful  —  namely,  the  study  of  the 
boundary-layer  properties  of  a  low  density  N2  nozzle  expansion  to  which  Na  is  added  in  very 
low  concent  rations  as  the  tracer  gas.  The  feasibility  of  applying  this  technique  to  such  a 
problem  was  investigated  by  developing  a  theoretical  description  of  the  time-resolved 
fluorescence  of  Na  in  a  N2  host  gas  environment,  and  the  results  of  this  work  are  described 
in  this  report. 

2.0  REVIEW  OF  LASER-INDUCED  FLUORESCENCE  (LIF)  MODELS 

The  basic  method  of  measuring  gas  species  parameters  by  the  laser-induced  fluorescence 
(LIF)  method  is  to  excite  an  atom  or  molecule  From  the  ground  electronic  state  to  an  excited 
electronic  state  by  absorption  of  visible  or  ultraviolet  laser  radiation.  The  subsequent 
fluorescence  produced  by  spontaneous  emission  back  to  the  ground  electronic  state  is  then 
measured.  The  intensity  of  the  fluorescence  is  proportional  to  the  total  number  density  of 
atoms  (or  molecules)  in  the  upper  excited  state  from  which  the  radiation  originates.  The 
total  number  density  in  the  original  excited  ground  state  is  then  related  to  the  upper  state 
number  density  by  an  appropriate  excitation-emission  model.  Finally,  the  total  species 
number  density  is  determined  from  an  appropriate  original  distribution  function  over  all  the 
possible  energy  states  which  may  be  populated  at  the  time  of  laser  excitation.  For  molecular 
species,  rhe  temperature  can  be  obtained  by  measuring  the  relative  fluorescence  intensities 
resulting  from  excitation  of  two  or  more  ground  slate  roiaLional-vibrational  levels.  For 
either  temperature  or  number  density  measurer  bnts,  the  important  step  in  obtaining  values 
of  those  parameters  w'ith  LIF  is  the  excitation-emission  model  used  to  relate  the  measured 
fluorescent  intensities  to  the  original  state  of  the  gas  species  being  investigated.  The  next  two 
sections  present  the  two  most  popular  models  used  in  LIF  analysis,  the  two-level  and  three- 
level  models. 

2.1  TWO-LEVEL  MODEL 

The  simplest  excitation-emission  model  is  the  two-level  model  shown  in  Fig.  I.  In  this 
model  it  is  assumed  that  the  only  important  energy  exchange  processes  are  those  which 
transfer  energy  between  the  two  laser-coupled  slates.  The  rate  equation  for  the  population 
of  level  2  can  be  written: 
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ill 


NIB|2'\  -  [Q21  '  A2I  I"  1*2,  V]  N2 


The  notation  for  this  and  subsequent  equations  is  described  in  the  Nomenclature. 


(I) 


Figure  1.  Two-level  model  schematic. 

If  it  is  assumed  that  at  time  =  0  all  of  the  population  is  in  level  I,  then 

nt  =  n,  +  n2  =  (2) 

For  steady-stale  conditions  the  rate  of  change  of  population  in  level  2  is  zero,  and  Eq.  (1) 
can  then  be  rearranged: 
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n 

1 


:  1 


(U 


t  2 


■*12^ 


(3) 


which  gives  the  ground  state  number  density  in  terms  of  known  (or  measurable)  constants 
and  the  number  density  in  level  2,  which  is  proportional  to  the  fluorescent  signal. 


If  the  laser  radiation  density  (e„)  is  large  so  that  (B]2  +  B21)  Qr  ^  Q21  +  A2i>  saturation 
conditions  exist,  and  Eq.  (2)  reduces  to 


which  allows  determination  of  the  gas  number  density  without  the  necessity  of  knowing  the 
quenching  rale  Q2|. 


Taking  Eq.  (2)  and  solving  for  N2  and  recognizing  that  the  fluorescent  signal  (S)  is 
related  to  the  upper  state  number  density  by 

S  .  C  A2I  Nj  (5) 
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then 


S  =  C  A 


KT  B12 


21 


^2]  '' 


A2i>/p„  +  <Bl2  +  D2|) 


(6) 


If  a  binomial  expansion  about  l/e„  is  performed,  the  following  expression  is  obtained: 


S  ,  C  A 


12 


21 


B12  "  H2  1 


N1  “ 


(7) 


A  plot  of  S  versus  i/p„  will  give  Nf  as  a  function  of  the  intercept,  and  0?i  as  a  function  of 
the  slope  (tor  known  N|).  This  “near”  saturation  technique  provides  another  method  for 
determining  the  number  density  without  knowledge  of  Q2|.  The  two-level  near-saturation 
technique  has  been  used  to  measure  number  densities  of  C2  (Ref.  1),  MgO  (Ref.  2),  CH  and 
CN  (Ref.  3),  and  Na  (Ref.  2)  in  flames.  Additionally,  OH  number  densities  have  been 
measured  using  the  two-level  model  at  low-pressure  short  pulse  conditions  to  eliminate 
quenching  effects  (Ref.  4).  Finally,  Bechtal  (Ref.  5)  has  used  the  iwo-level  model  with 
quenching  effects  included  for  measuring  OH  concentrations  in  a  flame. 


2.2  THREE-LEVEL  MODEL 


The  three-level  model,  as  shown  schematically  in  Fig.  2,  is  particularly  useful  when  two 
closely  lying  electronic  stales  occur  as  the  first  excited  levels  above  the  ground  state,  as  in 
sodium  and  potassium.  For  the  case  when  level  2  is  excited,  the  rate  equation  is 

dN2  r 

IT  '  ~LB21^+  A2J  '■  Q2I  H  Q23]N2H  [B|/„  »■  Qt2]Nl  +  Q32N,  (8) 


Assuming  that  all  of  the  population  is  in  level  1  before  laser  excitation,  then 

Nt  --  N,  +  N2  *  N,  =  NJ 
At  steady-stale,  Eq.  (6)  can  be  rearranged: 

N-n  =  _[(,I21  +  HQ12  .  y23]x2  -fg 32  -Q12-b]2PJn3 

B12PV  +  Q1  2 


(9) 


(10) 


If  all  quenching  and  collisional  exchange  rates  are  known,  then  Nj°  can  be  determined  from 
the  total  fluorescent  signal  from  levels  1  and  2.  Under  saturated  conditions  Eq.  (10)  reduces 
to 
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D 
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+ 


l  B.2  +  °2. 

\  B,2 


(ID 


which  would  allow  calculation  of  N°  from  LIF  measurements  withoul  the  necessity  of 
knowing  the  various  collisional  quenching  and  exchange  rates.  The  three-level  saturated  LIF 
technique  has  been  used  to  measure  Na  number  densities  in  a  flame  (Ref.  6). 


Figure  2.  Three-level  model  schematic. 

3.0  SODIUM-DIATOMIC  MOLECULE  SYSTEM  RELAXATION 


3.1  DETAILS  OF  QUENCHING  MECHANISM 


Experimental  (Refs.  7  through  9)  and  theoretical  (Ref.  10)  analysis  has  shown  that  the 
collisional  energy  exchange  collisions  between  sodium  and  diatomic  molecules  involve 
energy  transfer  between  the  electronic  states  of  sodium  and  the  vibrational  slate  of  the 
molecule.  The  quenching  (Qji,  Q31)  and  reverse  quenching  (Q12,  Q13)  rates  in  the  rate 
equations  for  the  sodium  levels  [Eq.  (8)]  are  therefore  dependent  on  the  vibrational 
population  distribution  of  the  quenching  molecule.  Specifically, 


On  -  £  ",  £ 

V  W  V 

On  -  £  Nw  E  K, 

W  V  <  vv 


(12) 


where  kvw  is  the  rate  constant  (cm3  sec-1)  for  an  upward  vibrational  transition  causing  a 
downward  sodium  transition  (quenching),  and  k'v  is  the  rate  constant  for  a  downward 
vibrational  transition  causing  an  upward  sodium  transition  (reverse  quenching).  The  upward 
and  downward  rate  constants  are  related  by  the  principle  of  detailed  balancing: 
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kwv  -  kvW^«p[-(AE2,-AEvJ/kT]  (13) 

where  g2  and  g,  and  the  statistical  weight  of  sodium  levels  2  and  1  (Fig.  2),  AE2,  is  the  sodium 
level  energy  difference,  and  AEV*  is  the  vibrational  transition  energy. 

During  an  excitation-relaxation  process  in  a  sodium-diatomic  mixture,  the  population  in 
the  vibrational  levels  of  the  quenching  molecule  may  not  maintain  a  Boltzmann  distribution 
because  of  the  relatively  slow  V-V  exchange  rates  (Appendix  A)  and  unequal  sodium- 
molecule  rate  constants  kvw  for  different  vibrational  transitions.  Consequently,  the  values  of 
the  total  quenching  and  reverse  quenching  rates  in  Eqs.  (8)  and  (10)  may  vary  with  time. 

3.2  SODIUM-NITROGEN  RELAXATION  MODEL 

To  check  the  effects  of  the  possible  variation  in  quenching  rales,  a  computer  code  was 
developed  to  calculate  the  time-dependent  populations  in  all  sodium  electronic  and 
molecular  vibrational  levels  during  laser  excitation  of  one  of  the  sodium  3S-3P  transitions 
for  nitrogen-sodium  mixtures. 


The  rate  equations  for  the  ground  electronic  state  (32S|/2)  and  doublet  first  excited  slates 
(32P1/2  and  32P2/2)  of  sodium  are,  for  excitation  of  the  32Pi/2  state  (level  2): 
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where  k?2  and  k23  are  the  P-level  exchange  rale  constants  for  collisions  with  sodium  atoms  in 
the  ground  state  (32Si/2-level  1)  and  k(2  and  k23  are  the  P-level  exchange  rate  constants  for 
collisions  with  nitrogen  molecules  (total  nitrogen  number  density  =  Nj).  It  is  assumed  that 
the  quenching  rate  constants  (k*v  and  are  equal  for  both  P-levels,  and  the  reverse 
quenching  rate  constants  from  detailed  balancing  are  calculated  [(Eq.  (13)]  so  that  for  v  <  w 

kiv  ^  kw> . 

For  a  nitrogen  vibrational  level  the  rate  equation  broken  into  its  component  parts  is 
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The  V-V  and  V-T  rate  constants  are  discussed  in  Appendix  A. 

If  the  complete  system  of  sodium  electronic  level  and  nitrogen  vibrational  level  rate 
equations  and  initial  conditions  of  a  Boltzmann  distribution  in  the  nitrogen  vibrational 
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levels  and  all  sodium  population  in  the  ground  electronic  state  is  used,  the  time-dependent 
populations  are  found  by  numerically  solving  the  system  of  equations  by  the  Gear  linear 
multistep  method  (Ref.  LI). 

Table  1  lists  the  values  of  the  energy  exchange  cross  sections,  along  with  some  of  the 
other  atomic  and  molecular  constants  used  in  this  analysis.  The  rate  constants  for  use  in  the 
rate  equations  arc  calculated  by  multiplying  the  cross  section  by  the  relative  speed  of  the 
colliding  species,  which  is  dependent  on  the  gas  temperature.  The  quenching  rate  constants 
are  based  on  the  change  in  vibrational  level  (Av)  during  the  quenching  collision.  A  Gaussian- 
like  distribution  of  the  quenching  probability  is  used,  as  suggested  by  the  experimental 
scattering  experiments  of  Ref.  12.  Figure  3  shows  the  probability  distribution  used  in  the 
present  work.  Notice  that  the  peak  is  placed  at  v  =  3  and  v  =  4,  not  at  the  resonance 
condition  (~v  =  8).  It  is  assumed  that  there  is  no  dependence  on  the  initial  (or  final) 
vibrational  level.  The  actual  quenching  rate  constants  are  calculated  by  taking  Ihe  relative 
probabilities  and,  using  a  Boltzmann  distribution  at  the  gas  temperature,  summing  the 
products  of  number  density  and  quenching  probability  for  each  vibrational  level  and  then 
normalizing  so  thaL 


£  \  £  i,„  -  q-i  cs) 

V  w  >  v 


Table  1.  Constants  for  Sodium-Nitrogen  Model 


— - 

Parameter 

Value 

Source 

A21  A31 

6.13  x  107  sec-1 

Ref.  17 

q32 

23  3  A2 

Ref .  1 8 

q23 

552  A2 

Detailed  Balance 

*  * 
q23 

144  A2 

Ref .  1 4 

>  * 
q32 

73.3  A2 

Detailed  Balance 

A(32P1/2  -  32S1/2) 

5895.9  A2 

Ref.  19 

'f3  ?3/2  "  3  S1/21 

5889.9  A 

Ref.  19 

*at  40 0 ° K 
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Figure  3.  Sodium-nitrogen  collisional  quenching 
probability  distribution. 

where  Qr  is  the  total  quenching  rate  (sec-  '),  This  rate  is  calculated,  in  turn,  from  a  chosen 
total  quenching  cross  section  and  the  relative  speed  of  the  colliding  species.  The  quenching 
rale  constants  for  the  sodium  32P.y2  and  3:Pj/2  levels  are  assumed  to  be  equivalent.  The 
reverse  quenching  rate  constants  are  calculated  from  the  quenching  rate  constants  by 
detailed  balancing  Eq.  (13).  The  reverse  quenching-rate  constants  for  the  two  sodium  P 
levels  will  be  similar  but  somewhat  different  in  value  because  of  the  17.3-cm-1  energy 
difference  between  the  two  P  levels,  and  to  the  different  statistical  weights  of  the  levels.  The 
sodium-nitrogen  total  quenching  cross  sections  (q_n  =  Q21 )  are  taken  from  the  review  article 
by  Lijnse  and  Elscnaar  (Ref.  8),  which  shows  an  inverse  temperature  dependence  of  the  total 
quenching  cross  section  below  about  1,000°K. 

Figure  4  presents  the  rate  constants  for  quenching  and  reverse  quenching  collisions 
between  the  sodium  3P  levels  and  nitrogen  for  vibrational  transitions  originating  (for 
quenching  collisions)  in  the  Olh  vibrational  level,  at  a  temperature  of  2,000°K.  From  the 
figure  it  is  apparent  that  reverse  quenching  rates  are  small  compared  to  quenching  rates  lor 
the  lower  v  levels  but  that  at  the  higher  v  levels  reverse  quenching  collision  rates  arc- 
comparable  to  quenching  rates.  As  shall  be  seen  in  the  next  section,  this  relationship  between 
the  quenching  and  reverse  quenching  rates  for  different  vibrational  levels  has  a  strong 
influence  on  the  time  history  of  the  laser-excited  sodium-nitrogen  system. 


Rate  Constant  (cm5  sec 
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Figure  4.  Quenching  and  reverse  quenching  rate  constants 
at  T  =  2,000°K. 
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The  3P  level  mixing  exchange  cross  sections  Tor  collisions  with  molecules  (qii,  qjj)  are 
temperature  adjusted  using  an  exponential  lit  to  the  velocity-dependent  cross  section  data  of 
Ref.  13  along  with  the  cross  section  at  4003K'from  Ref.  14.  Figure  5  show.*  the  cross  section* 
and  rate  constants  resulting  from  this  procedure.  Also  shown  in  the  figure  is  the  rate 
constants  which  would  result  if  a  constant  3P  mixing  cross  section  were  used.  The 
application  of  a  temperature  correction  on  the  3P  mixing  cross  section  reduces  the 
temperature  dependence  of  the  rate  constant  so  that  the  rate  constants  at  different 
temperature  are  more  nearly  equal.  The  3P  level  mixing  exchange  cross  sections  for 
collisions  with  3S  sodium  atoms  (q23,  q.32)  are  not  temperature-corrected  because  of  a  lack  of 
information  on  the  velocity  dependence  of  these  cross  sections.  However,  the  contribution 
of  the  3S  level  collisional  3P  mixing  exchange  rates  on  the  relaxation  in  a  sodium-nitrogen 
mixture  is  small  because  of  the  much  lower  sodium  numbei  densities  compared  to  the 
nitrogen  concentration  for  the  conditions  of  imciest. 


Figure  5.  Sodium  3P  level  mixing  exchange  rate  constants  and 
total  cross  sections. 

3.3  RESULTS  OF  SODIUM-NITROGEN  MODELING 

For  relatively  high  number  densities  of  nitrogen  and  sodium  (2  x  10|g  cm-*  and 
1.0  x  lO'-’-cm--1.  respectively)  and  a  temperature  of  l,000°K,  the  lime  history  of  the 
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populations  in  the  three  sodium  levels  and  in  the  nitrogen  vibrational  levels  is  shown  in  Figs. 

6  and  7  for  a  constant  laser  radiation  density  (step  function  on  at  t  =  0)  of  1.0  x  I0~l0ergs 

cm-3  Hz  - 1  exciting  the  3P1/2  level.  At  this  laser  excitation  energy,  which  is  well  below  the 
saturation  energy,  the  sodium  level  populations  reach  an  equilibrium  value  after  about  0.1 
nsec,  and  maintain  this  constant  level  until  approximately  0.1  msec,  when  a  transition  to  a 
new  steady-state  level  occurs.  The  vibrational  population  distribution  continually  changes  as 
more  and  more  molecules  arc  put  into  higher  vibrational  levels,  because  of  quenching 
collisions  with  sodium.  The  relatively  slow  V-V  rates  prevent  the  vibrational  state  from 
maintaining  an  equilibrium  distribution,  and  eventually  there  is  sufficient  change  from  the 
original  Boltzmann  distribution  to  alter  the  total  quenching  cross  section  and  to  permit 
significant  contributions  from  reverse  quenching  processes.  For  lower  temperatures  and/or 
higher  laser  intensities  the  results  at  these  high  number  densities  are  qualitatively  the  same. 
For  higher  temperatures  this  is  not  the  case,  as  illustrated  in  Figs.  8  and  9  for  identical 
conditions  to  those  in  Figs.  6  and  7  except  for  a  lemperature  of  2,000°K.  At  this  temperature 
there  is  no  transfer  to  a  second  steady-state  condition.  Despite  the  fact  that  there  is  a 
significant  change  in  the  vibrational  population  distribution,  this  change  is  not  large  enough 
to  cause  the  transfer  to  a  second  steady-state  condition.  The  vibrational  population 
distribution  attains  a  steady  state  after  approximately  1  msec,  so  that  there  is  no  possibility 
ot  a  change  at  times  longer  than  those  for  which  the  calculation  was  performed. 


Figure  6.  Time  history  of  population  in  sodium  levels  for 
high  density,  T  =  1 ,0Q0°K. 


16 


AEDC-TR-80-47 


Nitrogen 
Run  Number  58 

Mol.  N.D.  =  2.000  x  1019  cm-3 
Sodium  N.D.  ”  1.000  x  1013  m"3 
Temp.  -  1,000°K 
Num.  Vib .  Levels  =11 
3P,  Excited 

Laser  Rad.  Don.  =  1.0  x  10 
Total  Q  “  20  Sq.  Ang. 


Figure  7.  Time  history  of  population  in  nitrogen  vibrational  levels 
for  high  density,  T  =  1 ,000°K. 


At  lower  number  densities  a  still  different  picture  emerges.  Figure  10  shows  the  time 
history  results  for  nitrogen  and  sodium  number  densities  of  2.0  x  10i7  and  1  x  I010  cm--,  a 
temperature  or  2,000°K  and  laser  radiation  density  of  1.0  x  10-10  ergs  cm-3  Hz-1.  The 
steady-stale  condition  is  not  attained  until  after  about  7  nsec  with  the  3P|/2  level  showing  a 
pronounced  peak  at  about  0.2  nscc.  This  behavior  occurs  because  of  the  dependence  of  the 
3P  level  mixing  exchange  rates  on  the  nitrogen  number  density.  For  these  lower  number 
densities  the  3P3/2  populating  collisions  occur  al  a  reduced  rate  so  that  the  3P3/2  level 
population  process  will  lag  the  laser-excited  3P|/2  level  population.  Al  these  conditions  the 
3P|/2  level  attains  a  number  density  higher  than  its  steady-stare  value  before  collisions  with 
nitrogen,  causing  a  transfer  lo  the  3P3/2  level  to  become  effective.  After  about  0.2  nsec  this 
transfer  process  begins  to  depopulate  the  3P1/2  level  until  the  sready-state  condition  is 
obtained.  At  lower  temperatures  (Fig.  1 1)  the  same  picture  emerges  with  the  exception  that 
at  very  long  times  (-0.1  sec)  a  transition  to  a  second  steady-state  condition  occurs,  because 
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Figure  8.  Time  history  of  population  in  sodium  levels  for 
high  density,  T  =  2,000°K. 

of  the  redistribution  of  population  in  the  vibrational  levels  of  nitrogen  (Fig.  12).  It  takes  a 
longer  time  at  these  conditions  to  attain  sufficient  population  in  the  higher  vibrational  levels 
to  cause  the  transfer  to  a  second  steady  state  than  for  the  higher  number  density  conditions. 
It  should  also  be  noted  that  at  the  same  laser  radiation  density  the  system  is  closer  to 
saturation  (N|  =  N2)  at  the  lower  nitrogen  densities  because  of  a  reduction  in  the  total 
quenching  rale. 

In  summary,  this  analysis  has  shown  that  there  can,  in  fact,  be  a  significant  effect  on  the 
relaxation  of  a  sodium-nitrogen  laser-excited  system  attributable  to  redistribution  of 
vibrational  population  by  quenching  collisions.  However,  this  effect  occurs  at  sufficiently 
long  times  so  that,  for  short  laser  pulse  measurements  (<1  gscc)  the  steady-slate  relations 
[Eq.  (10)]  for  a  three-level  system  arc' applicable.  On  the  other  hand,  extremely  short  pulses 
may  not  be  long  enough  for  steady  state  to  be  obtained  under  some  conditions  because  of 
the  finite  3P  level  exchange  collision  rate. 
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Figure  9.  Time  history  of  population  in  nitrogen  vibrational  levels 
for  high  density,  T  =  2,000°K. 
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Figure  11.  Time  history  of  population  in  sodium  levels  for 
low  density,  T  =  300° K. 
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Figure  12.  Time  history  of  population  in  nitrogen  vibrational  levels 
for  low  density,  T  =  300°K. 


AE0C-TR-80-47 


A  E  DC-TR  -80-47 


4.0  SINGLE  DIATOMIC  SPECIES  RELAXATION 

4.1  BREAKDOWN  OF  TWO-LEVEL  MODEL  ASSUMPTIONS 

Real  molecules  contain  complex  rotational  and  vibrational  structure  which  may  cause 
significant  breakdown  of  the  assumptions  in  the  two-level  L1F  model.  Collisional  energy 
exchange  collisions  during  the  laser  excitation  pulse  will  transfer  energy  between  the  laser- 
coupled  vibrational-rotational  states  and  the  other  vibrational-rotational  states  of  the  upper 
and  lower  electronic  states  at  a  finite  rate.  Additionally,  spontaneous  emission  from  the 
laser-populated  upper  level  will  occur  not  only  to  the  laser-excited  lower  level,  but  also  to 
other  vibrational  levels  of  the  ground  electronic  state.  Quenching  collisions  will  also  produce 
transitions  to  lower  electronic  vibraiional-rotational  states  other  than  the  laser-excited  one. 
The  possible  effects  of  Lhcse  processes  will  be: 

1.  Time-dependent  populations  and  fluorescent  intensities 

2.  Non-Boltzmann  vibrational  and  rotational  population  distributions 

3.  Time  dependent  rotal  number  density  in  laser-coupled  states  [breakdown  of  Eq. 

(2)]. 

Luchl  and  Laurendeau  (Ref.  15)  have  shown  that  when  the  ratio  of  the  rotaiional  relaxation 
rates  of  the  upper  and  lower  electronic  states  to  the  quenching  rate  are  100  or  less,  the 
steady-state  rotational  state  population  distributions  are  significantly  non-Boltzmann.  To 
investigate  the  steady-state  vibrational  population  distribution  and  the  time  history  of  the 
rotational  and  vibrational  distributions,  a  computer  code  has  been  constructed  which 
follows  the  detailed  relaxation  of  a  diatomic  molecule  during  laser  excitation.  The  next 
section  discusses  the  model  used  in  the  computer  code. 

4.2  SINGLE  DIATOMIC  SPECIES  RELAXATION  MODEL 

For  the  case  of  a  single  diatomic  species  in  which  one  particular  vibrational-rotational 
( V-R)  level  in  the  ground  electronic  state  of  the  molecule  is  pumped  by  laser  radiation  to  one 
particular  vibrational-rotational  level  in  the  first  excited  electronic  stale  of  the  molecule,  the 
following  processes  need  to  be  considered: 

1.  Laser  excitation  and  stimulated  emission 

2.  Fluorescence  from  all  upper  electronic  V-R  states 
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3.  Collisional  exchange  between  rotation  and  translation  (R-T) 

4.  Collisional  exchange  between  the  vibrational  modes  of  collision  partners  (V-V) 

5.  Collisional  exchange  between  vibration  and  translation  (V-T) 

6.  Collisional  quenching,  i.e.,  nonradiative  transitions  from  upper  to  lower 
electronic  state 


7.  Radiative  transitions  between  V-R  levels  of  the  same  electronic  state 
(heteronuclear  molecules  only) 

8.  Self-absorption  of  fluorescence 


It  is  assumed  that  the  gas  is  optically  thin  so  that  process  8  can  be  neglected.  Additionally, 
the  radiative  lifetime  for  transitions  between  V-R  levels  of  the  same  electronic  state  is 
typically  much  greater  than  the  time  scales  for  processes  1  through  6,  so  process  7  will  also  be 
neglected.  Including  processes  1  through  6,  the  rate  equation  for  the  lower  electronic  level, 
laseT-cxcited  V-R  state  is  then  (with  its  component  parts  specified): 
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populating  V  I’  &  R  I’ 


(17) 


Where  v"  and  J"  denote  the  laser-excited  lower  electronic  state  vibrational  and  rotational 
level,  v*  and  J'  denote  the  laser-excited  upper  electronic  state  vibrational  and  rotational 
level,  N|  and  N2  are  the  lower  and  upper  electronic  slate  number  densities  respectively,  NT  is 
the  total  gas  number  density,  B|2  and  B21  are  the  absorption  and  simulated  emission 
coefficients,  e,,  is  the  laser  radiation  energy  density,  A  the  radiation  rate  (Einstein 
coefficient),  Q  the  quenching  rate,  and  the  R’s  are  the  various  collisional  exchange  rate 
constants  (cm-Vsec)  which  are  identified  by  their  subscripts.  Rvw  is  the  V-T  rale  constant 
for  transfer  from  vibrational  state  v  to  vibrational  state  w'.  Ru;  is  the  R-T  rate  constant  for 
for  transfer  from  rotational  state  J  to  state  J',  and  Rv+)  «,■.  v,  w  +  i  is  the  rale  constant  for 
V-V  exchange  for  molecules  in  vibrational  states  v+1  and  w'  colliding  and  exchanging 
vibrational  energy  to  end  up  in  states  v  and  w+  1.  For  V-T  and  V-V  processes,  |Av|  =  1 
transitions  only  are  considered,  but  for  R-T  processes  all  AJ’s  are  included  in  the  analysis. 
The  sum  over  J  for  the  A  and  Q  terms  is  limited  to  allowed  dipple  transitions  only,  i.e.,  P, 
Q,  and  R  branch  transitions  in  the  case  of  CO.  The  quenching  rates  are  based  on  the  value  of 
Lhe  Einstein  A  coefficients  for  the  various  vibrational-rotational  transitions. 


The  rate  equation  for  the  laser-excited  upper  electronic  V-R  state  is 
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(18) 


The  rate  equations  for  the  other  V-R  levels  in  the  upper  and  lower  electronic  states  are  the 
same  as  those  for  the  laser-coupled  states  without  the  absorption  and  stimulated  emission 
terms. 
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4.3  CARBON  MONOXIDE  MODELING  RESULTS 

The  total  number  of  vibrational-rotational  levels  used  in  the  relaxation  model  was 
restricted  by  computer  time  and  core  usage  limitations;  it  was  not  possible  to  use  a  large 
number  of  rotational  and  vibrational  levels  simultaneously.  The  analysis  was  therefore 
divided  into  two  parts.  The  first  part  used  31  rotational  levels  and  3  vibrational  levels  to 
investigate  the  rotational  relaxation,  and  the  second  part  neglected  rotational  relaxation  and 
used  10  vibrational  levels  in  each  electronic  state.  For  both  test  cases  a  step  function  laser 
pulse  was  applied  at  time  t  =  0  and  a  constant  laser  radiation  density  was  maintained  for  the 
duration  of  the  run  time.  Excitation  is  from  a  rotational  level  in  the  0th  vibrational  level  of 
the  X  '£+  ground  electronic  state  of  CO  up  to  one  particular  vibrational-rotational  stale  in 
the  A'n  first  excited  electronic  state.  The  absolute  transition  probabilities  (Einstein  A 
values)  of  Ref.  16  are  used  for  the  various  vibrational  transitions  from  the  upper  to  the  lower 
electronic  state. 

The  time  history  of  the  rotational  population  distributions  in  the  upper  and  lower  laser- 
coupled  vibrational  levels  is  shown  in  Figs.  13  and  14  for  atmospheric  pressure  conditions  at 
300°K  and  a  laser  radiation  density  of  1 .0  x  10-8  ergs  cm -3  Hz~  >,  which  is  well  under  the 
saturation  radiation  density.  The  results  agree  qualitatively  with  those  in  Ref.  17.  The 
distributions  are  highly  non-Boltzmann  with  an  overpopulation  in  the  upper  electronic  state 
laser-pumped  rotational  level  and  an  underpopulation  in  the  lower  electronic  state  laser- 
pumped  rotational  level. 


Rolutinna1  Qjunfum  Number. .  J 

Figure  13.  Time  history  of  rotational  population  distribution  in  lower 
electronic  state  laser-coupled  vibrational  level. 
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Carbon  Monoxide 

Total  Gas  N.D.  -  2.448  x  1019  cm'9 
Gas  Temp.  -  300°K 
Nun.  Rot.  Levels  -  31 
Num.  Vib.  Levels  =  3 

J  =  6  V  -  0  to  J  ■  5  V  =  1  Excitation 
Laser  Bad.  Den.  *  1.0  x  10-8 
Q  -  100  x  A 
Lower  Elec.  Lovol 


Rotational  Quantum  Number,  J 


Figure  14.  Time  history  of  rotational  population  distribution  in  upper 
electronic  state  laser-coupled  vibrational  level. 


For  pure  vibrational  relaxation,  Fig.  15  shows  the  time  history  of  the  number  density  of 
the  upper  electronic  slate  laser-excited  vibrational  level  for  room  temperature,  atmospheric 
pressure  conditions  with  three  values  of  the  quenching  rate  (Q).  For  all  cases  the  population 
does  not  come  to  steady  state  until  very  long  times  (—0.1  sec).  The  initial  drop  in  number 
density  (around  0.1  nsec  for  Q  =  100  x  A)  depends  strongly  on  the  quenching  rate  and  is 
attributable  to  the  depletion  of  population  from  the  upper  laser -excited  state  by  quenching 
and  radiative  relaxation  to  the  various  ground  electronic  slate  vibrational  levels  where  the 
finite  V-V  rates  prevent  fast  enough  collisional  exchange  to  refill  the  depleted  laser-excited 
vibrational  level.  Figures  16  and  17  show  the  lime  history  of  the  vibrational  population 
distributions  in  the  upper  and  lower  electronic  states  for  Q  =  100  x  A.  Highly  non- 
Boltzmann  distributions  exist  at  all  times  during  the  relaxation.  V-V  collisional  exchange 
tends  to  smooth  out  the  vibrational  population  distribution  in  the  ground  electronic  state, 
with  the  most  pronounced  effects  occurring  after  about  1  msec.  At  a  temperature  of 
2,000° K,  (Fig.  18)  the  time  history  is  qualitatively  similar  to  that  at  300°K  except  for  a  higher 
absolute  steady-state  number  density  at  an  earlier  time.  This  effect  is  attributable  to  the 
much  higher  V-V  and  V-T  rates  at  the  higher  temperature. 
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Figure  15.  Time  history  of  population  in  upper  electronic  state  laser-coupled 
vibrational  level  for  three  values  of  quenching  rate. 
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Figure  16.  Time  history  of  vibrational  population  distribution  in 
upper  electronic  state,  T  =  300“K. 
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Figure  17.  Time  history  of  vibrational  population  distribution  in 
lower  electronic  state,  T  =  300°K. 
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Figure  18.  Time  history  of  population  in  upper  electronic  state 
laser-coupled  vibrational  level,  T  -  2,000°K. 
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The  effects  of  laser  radiation  density  on  the  upper  vibrational  level  number  density  arc 
presented  in  Fig.  19  and  Table  2.  In  the  figure  is  plotted  the  peak  number  density  in  the 
upper  vibrational  level  normalized  by  the  steady-state  number  density  predicted  by  the  two- 
level  model  for  both  300°K  and  2,000°K.  At  the  lower  temperature  the  peak  number  density 
is  equal  to  the  two-level  model  prediction  for  laser  radiation  densities  greater  than  10-6  erg 
env1  Hz'1.  The  ratio  falls  quickly  around  I  O  ’*  erg  cm  -  H/-'  and  again  levels  out  at  a  lower 
lower  ratio  for  radiation  densities  less  than  I0_]0  erg  cm  -1  Hz-1.  For  the  higher 
temperature  case  the  same  pattern  emerges,  except  that  at  the  saturated  conditions  the  ratio 
is  below  a  value  of  one.  In  Table  2  the  steady-state  time  and  number  density  are  essentially 
constant  for  all  the  laser  radiation  densities  tested  at  300°K.  At  2,000&K  the  time  of  steady 
state  is  virtually  constant  but  the  number  density  does  show  a  reduced  value  at  the  lower 
radiation  densities. 

In  summary,  for  saturated  laser  excitation  of  carbon  monoxide  the  peak  of  the  number 
density  in  the  upper  laser-excited  vibrational  level  is  comparable  to  the  two-lcvcl  prediction 
at  room  temperature.  For  lower  radiation  densities,  for  high  temperatures,  or  for  prediction 
of  the  lime  history  and  steady-state  conditions,  the  iwo-level  model  fails  completely  because 
oT  the  finite  collisional  exchange  rates  and  large  contribution  of  multipath  radiative  and 
collisional  relaxation  from  the  upper  electronic  stale  vibrational  manifold  to  the  ground 
electronic  state  vibrational  levels. 


Figure  19.  Laser  radiation  density  effects  on  ratio  of  calculated  peak 

number  density  and  two-level  model  steady-state  prediction. 
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Table  2.  Maximum  and  Steady-State  Number  Densities  and  Times 
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5.0  SUMMARY 

The  objective  of  this  work  was  10  assess  ihc  feasibility  of  using  the  laser-induced 
fluorescence  technique  to  measure  space-and  lime-resolved  species  number  densities  in 
gaseous  systems.  Specifically,  computer  codes  were  developed  to  model  the  detailed  kinetic 
rate  relaxation  of  Na  atom — N2  molecule  and  CO  molecule  systems  to  evaluate  the 
usefulness  ol  the  simple  two-  and  three-energy  level  models  in  deriving  number  densities 
from  fluorescence  intensities. 

For  sodium-nitrogen  mixtures  the  time-dependence  of  the  fluorescence  intensity  was 
shown  to  be  dependent  upon  both  the  Lotal  N2  number  density  and  the  N2  vibrational  mode 
number  density'  distribution  function.  Specifically,  it  was  shown  that  a  mctastable  steady- 
staLe  fluorescence  intensity  of  Na  was  attained  within  0.1  -  10  nsec  after  initiation  of  the 
laser  source  radiation.  Further,  the  metastable  steady-state  fluorescence  intensity  for  the 
three-level  Na  model  was  dependent  upon  the  N2  number  density,  thereby  providing  a  means 
of  measurement  of  the  local  N2  number  density  using  Na  as  a  fluorescent  tracer.  As  time 
increases,  the  effects  of  vibration-vibration  and  vibration-translation  energy  transfer 
processes  in  the  N2  assume  added  importance,  and  a  transition  occurs  in  the  fluorescence 
signal  from  the  metastable  state  to  a  final  steady  state.  For  the  parameters  used  in  this  study 
it  was  seen  that  this  transition  occurred  approximately  1  msec  after  initiation  of  the  laser 
source  radiation.  Consequently,  for  laser  pulses  sufficiently  long  in  duration  to  attain  the 
metastable  steady-state  condition,  and  shorter  than  the  time  of  transition  to  the  final  steady- 
state  condition,  the  laser-induced  fluorescence  technique  using  the  three-level  model  appears 
to  be  feasible  for  measuring  nitrogen,  or  other  similar  diatomic  molecule,  number  densities 
with  sodium  as  a  tracer.  Experimental  measurements  under  controlled  conditions  are 
recommended  to  validate  further  the  applicability  of  this  technique. 

For  molecular  systems  the  preliminary  results  using  pure  carbon  monoxide  have  shown 
the  vital  importance  of  multipath  relaxation  from  the  upper  laser-excitcd  level.  For  carbon 
monoxide  the  two-level  model  is  totally  inadequate  for  predicting  the  steady-slate 
populations.  Consequently,  the  laser-induced  fluorescence  technique  based  on  the  two-level 
model  cannot  be  used  with  any  degree  of  reliability  for  carbon  monoxide  number  density 
measurements.  Furthermore,  these  results  make  it  clear  that  it  is  highly  desirable  to  extend 
the  computer  modeling  techniques  developed  here  to  more  realistic  multispccies  molecular 
systems  and  to  species  of  interest,  such  as  OH,  in  order  to  evaluate  the  applicability  of  laser- 
induced  fluorescence  to  these  systems. 
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NOMENCLATURE 

AtJ  Einsiein  spontaneous  emission  coefficient  rate  for  i-j  transition 

B,j  Radiative  absorption  coefficient  (i  <  j),  stimulated  emission  coefficient 

(i  >  j) 

C  Constant  proportional  to  viewing  volume  and  optical  parameters 

[(Eq.  (5)] 

E  Energy 

gj  Statistical  weight  of  level  i 

k  Boltzmann’s  constant 

ktJ  Rate  constant  involving  energy  slates  i  and  j 

N,  Number  density  of  level  i 

N'r  Total  diatomic  ga.s  species  number  density 

Q,j  Quenching  (i  >  j)  or  reverse  quenching  (i  <  j)  rate 

Qt  Total  quenching  rate 

q,j  Cross  section  for  i-j  collisional  exchange 

Rl>);v.w  V-V  rate  constant  for  transition  from  stales  i,  j  to  stales  v,  w 

Rij  R-T  or  V-T  rate  constant  for  i-j  transition 

S  Fluorescent  signal 

T  Kinetic  temperature 

t  Time 

Vibrational-vibrational  transitions 


V-V 


Vibrational-rotational  transition 


Vibrational-translational  transition 
Vibrational  quantum  level 
Wavelength 


Excitation  radiation  energy  density 
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APPENDIX  A 

MOLECULAR  EXCHANGE  RATE  CONSTANTS 
1.0  INTRODUCTION 

The  theoretical  formulation  of  Herzfeld  and  Litovitz  (HL)  {Ref.  A-l)  for  collisiona! 
energy  exchange  attributable  to  short  range  interactions  is  employed  to  calculate  molecular 
V-V  and  V-T  exchange  rate  constants.  Their  expressions  yield  values  for  the  V-T  and  V-V 
transition  probabilities  (P)  based  on  the  kinetic  (translational  energy  exchange)  rate  constant 
(R()  at  a  given  temperature,  so  that  the  V-T  or  V-V  rate  constant  will  be 

II 

the  kinetic  collision  rate  (Rt)  is 


where  d  is  the  collisional  interaction  radius  (molecular  diameter),  c,  is  the  average  relative 
velocity  of  I  he  colliding  molecules,  a  the  collision  cross  section,  and  m  is  the  reduced  mass. 

2.0  VIBRATIONAL  EXCHANGE 


=  PR. 


M' 


(A- 1 ) 


(A-2) 


2.1  V-T  PROBABILITIES 


Following  the  HL  analysis,  the  V-T 
components 


transition  probability  is  separated  into  three 


P 


v  r 


JL  1  i 

Z_  /.  u' 


(A-3) 


Z,)  is  a  geometrical  factor  estimated  by  HL  to  be  3.  Single  quantum  vibrational  transitions 
only  are  considered  in  (he  HL  formulation.  For  (he  general  case  of  vibration  exchange  in 
molecule  B-C  attributable  to  a  V-T  collision  with  molecule  A, 


( MjJ  ¥  Mf“  )  M, 


Z,.vC 


My  Mj.  tMj 


Mu  r  Mc ) 


(A-4) 


where  M.\,  M|3,  and  Mc  are  the  molecular  weights  of  molecule  A,  atom  B,  and  atom  C 
respectively,  0  is  the  characteristic  temperature  of  the  vibrational  transition, 
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for  vibrational  energies  E|  and  E2,  and  finally 

a  ■  _  1  fi  ff4  m  I:  “  v* 

k 


(A-6) 


where  m  is  the  reduced  mass  of  the  system 


V1,v  mhi: 


ma  -  mhc 


(A-7) 


v  is  the  frequency  corresponding  to  the  vibrational  transition  (j<  =  Ej-E^/h),  k  is 
Boltzmann's  constant  and  f  is  the  range  of  repulsive  forces  between  the  colliding  molecules. 
Next 


] 


Zvr' 


for  translational  temperature  T. 


(A-8) 


The  final  form  of  the  HL  V-T  probability  is 


P 


VC  - 


J_  _l _ 1_ 

Zo  tost,  Zt.r 


Y  (2,  2)  c'f/kr(l 


(A-9) 


where  rc/r0  is  the  ratio  of  the  stopping  point  (rc)  for  a  molecule  coming  from  infinity  and 
undergoing  an  inelastic  collision  and  the  Lennard-Jones  constant  (r0),  Y(2,2)  is  the 
“Sutherland  correction”  factor 


Y  (2,  2)  =  0.76  ^1  r  1.1  (A-10) 

and  e  is  the  Lennard-Jones  potential  well  depth.  Table  A-l  lists  the  values  for  the  constants 
used  in  the  sodium-nitrogen  and  carbon  monoxide  relaxation  codes. 

2.2  V-V  PROBABILITIES 


The  V-V  probability  is  similarly  divided  into  three  components: 
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Table  A-1.  Constants  for  V-T,  V-V  and  R-T  Calculations 


where 


Parameter 

Value 

Source 

0 

3.75 

Ref .  A-4 

r/r 
c  o 

1.0 

Approximation 

e/k 

91.5 

Ref.  A-1 

£ 

2.117  x  10~9 

Ref.  A-1 

F 

0.15 

* 

C 

0.10 

* 

A 

See  Appendix  Section  3.0 
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(A- 12) 


with 


5,2  - 


16 jt4  m12  k  P2  (0,  -02)2 


(A- 1 3) 


Note  that  the  third  term  in  the  exponential,  because  of  the  anharmonicity  of  vibrational 
motion,  accounts  for  the  difference  between  endothermic  and  exothermic  collisional 
transfer,  i.e.,  detailed  balance  effects.  In  the  actual  calculation  of  1/Z,r  for  both  V-T  and 
V-V  probabilities,  exothermic  probabilities  only  are  calculated  using  the  full  equations,  with 
the  corresponding  reverse  probabilities  calculated  using  the  detailed  balance  relation 
Pfll  —  Pcx  g-AE/kT  where  AE  is  the  energy  given  up  to  translation  during  the  collision  because 
of  the  mismatch  in  vibrational  energies  during  the  V-V  exchange  collision.  Now, 
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Mn2  -  M®  \  6,6 


I  “2 

H  "'C  /  2  Q\  9 2 


4  A! M 


(A- 14) 


Substituting  for  0'  and  0  and  simplifying  yields 
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(A- 1 5) 
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Z0>  sometimes  referred  to  as  the  steric  factor,  is  the  term  which  adjusts  for  the  fact  that  the 
theoretical  analysis  which  produces  the  mathematical  expressions  for  Z,,  and  Z0,v  neglects 
the  rotational  motion  of  the  colliding  molecules.  HL  suggest  a  value  of  Z0  =  9. 

Figure  A- 1  shows  the  V-V  and  V-T  transition  probabilities  at  various  translational 
temperatures  for  nitrogen.  For  temperatures  below  approximately  1,000°K  the  V-T 
probabilities  for  the  lower  vibrational  levels  are  negligible  compared  to  the  V-V 
probabilities.  Despite  the  smallness  of  the  V-T  probabilities  it  was  necessary  to  include  V-T 
transfer  in  the  relaxation  code  to  obtain  steady-slate  conditions  at  a  finite  time. 

For  carbon  monoxide  it  is  known  that  long  range  dipole  interactions  play  a  significant 
role  in  V-V  exchange  processes  for  the  lower  vibrational  levels  (Refs.  A-3  and  A-4).  This 
effect  tends  to  increase  the  V-V  probabilities  calculated  by  the  HL  method  by  approximately 
two  orders  of  magnitude  at  a  temperature  of  300°K  for  collisions  involving  the  lower 
vibrational  levels  (v  <  10),  although  the  shape  of  the  probability  distribution  is  also  altered 
so  that  the  deviation  from  the  HL  probabilities  is  not  constant  for  each  transition. 
Nonetheless,  for  the  purposes  herein,  which  are  primarily  qualitative,  the  HL  V-V 
probabilities  multiplied  by  100  are  used  as  representative  of  the  true  V-V  probabilities  for 
carbon  monoxide. 


3.0  R-T  RATE  CONSTANTS 

The  exponential  gap  law  as  proposed  by  Polanyi  and  Woodall  (Ref.  A-2)  is  used  to 
calculate  the  R-T  probabilities.  For  a  downward  transition  from  J  to  J'  the  transition 
probability  is 

Pjj  -  =  F  exp  [^-C  (L  j  —  Ej')/k'fJ  (A-16) 

where  F  is  the  normalization  constant,  and  C  is  the  parameter  governing  the  exponential 
decrease  in  the  probability  with  increasing  energy  gap  AE  =  Ej  -  Ey. 

Table  A-l  gives  the  values  of  the  constants  in  Eq.  (A-16)  which  were  used  in  the  present 
analysis.  These  values  are  somewhat  arbitrary  because  of  the  lack  of  any  experimental  values 
for  the  R-T  exchange  rates  in  CO.  They  are  chosen  so  that  the  probability  distribution 
variation  with  AJ  and  J  approximates  that  found  for  HC1  in  Ref.  A-2.  The  value  of  Pjj-  for 
J  =  J '  (AJ  =  0)  is  adjusted  so  that  the  total  probability  is  equal  to  one. 
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NOMENCLATURE 

C  Exponential  parameter  [(Eq.  (A- 16)1 

cr  Average  relative  velocity 

d  Collision  interaction  radius  (molecular  diameter) 

Ej  Energy  of  level  i 

F  Normalization  constant  [Eq.  (A- 16)] 

h  Planck’s  constant 

k  Boltzmann’s  constant 

f  Repulsive  intermolecular  force  range 

M,  Molecular  weight  of  species  i 

in  Reduced  mass 

P  Transition  probability 
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R. 


rc 

ro 

T 

Y(2,2) 

Zo»Zosc*Z>ir 


e 

0 


(9.2' 


<T 


Collision  rate  defined  by  Eq.  (A-2) 

Stopping  point  in  molecular  collision 
Leonard- Jones  constant 
Temperature 

Sutherland  correction  factor 

Components  of  vibrational  transition  probabilities 

Lennard-Jones  potential  well  depth 

Characteristic  temperature 

Term  defined  by  Eq.  (A-6) 

Term  defined  by  Eq.  (A-13) 

Frequency 

Collisiona!  cross  section 
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